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hard vs. soft constraint 
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The curvature potential arising from confining a particle initially in three-dimensional space onto 
a curved surface is normally derived in the hard constraint q —* limit, with q the degree of freedom 
normal to the surface. In this work the hard constraint is relaxed, and eigenvalues and wave functions 
are numerically determined for a particle confined to a thin layer in the neighborhood of a toroidal 
surface. The hard constraint and finite layer (or soft constraint) quantities are comparable, but 
both differ markedly from those of the corresponding two dimensional system, indicating that the 
curvature potential continues to influence the dynamics when the particle is confined to a finite 
layer. This effect is potentially of consequence to the modelling of curved nanostructures. 

PACS numbers: 03.65Ge, 68.65.-k 



I. INTRODUCTION 

The existence of a potential Vc in the Schrodinger 
equation which stems from constraining a particle to a 
one or two-dimensional surface embedded in three dimen- 
sions has a long history 0, S S S H @ . The manifes- 
tations of Vc have been investigated through formal and 
numerical means 0, 13, S E3, Ell E3, E3, E3, EM [1^1 , 
motivated recently in part by the sophistication with 
which nanostructures can be fabricated. The physics 
of objects with novel geometries is increasingly relevant 
to the modelling of real devices, hence substantial effort 
has been directed towards understanding the physics of 
bent tubes and wires, as well as more complicated shapes 

ElilS El El El II. 

Consider a surface £(u,i>) with (u,v) surface coordi- 
nates and define q as the coordinate labelling the degree 
of freedom normal to E(u,w). Generally Vc (as detailed 
in the following section) is derived by imposing a hard 
constraint on the particle wherein a q — * limit is taken 
along with a wave function re-scaling such that the norm 
is preserved. Here, rather than imposing a hard con- 
straint, the particle will be confined to a thin layer in 
the neighborhood of a toroidal surface. The extent to 
which the hard constraint mirrors the more physically 
realizable soft constraint is then determined by calculat- 
ing some low-lying eigenvalues and eigenfunctions of the 
system. 

There are several reasons to investigate these ideas 
with a torodial structure: 

1. The symmetry of the torus reduces computational in- 
tensiveness, but, because a torus has non-trivial mean 
and Gaussian curvatures, curvature effects remain im- 
portant pn. 

2. The spectrum and eigenfunctions for a particle on a 
toroidal surface have been determined |2(j , so compar- 
isons can be made between the finite layer system and 
the two dimensional system both with and without Vc 
present. 

3. Toroidal structures have been fabricated and calcula- 



tions addressing their transport E3, EE E! E3 an d mag- 
netic [3l| properties performed. Toroidal structures are 
novel because unlike a bulk sample, conductivity through 
the device is anticipated to be dominated by azimuthal 
modes. 

The remainder of this paper is organized as follows: 
in section II, H qi the Hamiltonian for a particle near a 
toroidal surface is derived. The hard constraint q — > 
limit of that Hamiltonian is then taken; under the re- 
quirement that the norm of the wave function be pre- 
served, He obtains. Finally, the ab initio q = Hamil- 
tonian Hq is written. In section III the computational 
method employed to generate eigenvalues and wave func- 
tions is presented. Section IV gives results and section V 
is reserved for conclusions. 



II. THE TOROIDAL SCHRODINGER 
EQUATIONS 

To restate, there are three Hamiltonians relevant to 
this work: 

1. H q will be the Hamiltonian for a particle allowed to 
move in a thin layer normal to T 2 , where again the nor- 
mal degree of freedom will be labelled by q. 

2. He will be the Hamiltonian derived from H q after 
imposing the q — > hard constraint, and 

3. Hq the Hamiltonian for a particle restricted ab initio 
to T 2 , i.e., q = at the onset of the derivation of Ho- 
lt is best to begin with the most general case, that of 

H qi and later take the appropriate limits to obtain He 
and Hq. 

Points near a toroidal surface of major radius R and 
minor radius a may be parameterized in terms of cylin- 
drical coordinate unit vectors and a vector n normal to 
the surface by (2(| 

r(6, <p,q) = (R + a cosd)p + a sin6>k + qh. (1) 
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Applying d to Eq.(2) gives 

dr = (a + q)d99 +(R+(a + q) cosi 



dqn (2) 



with n = <p x 9, and 9 = —sm9p + cos#k. The metric 
elements gij can be read olf of 



dr ■ dr = (a + qfd9 2 + (R+(a + g)cos6 



dq 2 (3) 



and the Laplacian derived from 



_d_ 

dq 1 



9 2 9- 



(4) 



Setting a q = a + q and F q = R + (a + q)cos9 yields the 
H q Schrodinger equation 



1 d 2 ip sin# dip 1 d 2 ip 



a 2 d9 2 a q F q 89 



F| dcj) 2 



2h 



dip 



d 2 ip 



2V n (q) + 2EiJj = 



(5) 



r 



with h the mean curvature given by 



h= i(fci + fc 2 ) = \ 



1 



a„ 



cos9 



(6) 



It is convenient to also define the Gaussian curvature k 

1 cos# 



k 1 k 2 



TP 

q ?q 



(7) 



To derive the Schrodinger equation appropriate to He, 
V n (q) must be chosen to drive the particle arbitrarily 
close to the q = limit 0, H 0. As q — > the wave 
function is expected to decouple into surface and normal 
parts as 



^(9,(j ) ,q)^Xs(e,<j>)Xn(q). 



(8) 



Conservation of the norm must be preserved leading to 

Hi ii mm 



\iP\ 2 WdSdq= \ Xs \ 2 \Xn\ 2 dSdq 



1p = XsXnW ^ 



(9) 



(10) 



where W = 1 + 2qh + q 2 k and dS the surface measure. 
Performing the differentiations and letting q — > gives 
the pair of equations 



d 2 ip a sin 6 dip a 2 d 2 ip 



89 2 



F 89 F 2 



2a 2 V c + Pip = 0, (11) 



Making the standard ansatz for the azimuthal part of 
the eigenfunction x{4>) — ex P [im<j>] reduces Eq. (11) to 



d 2 ip 
~dff 2 



a sin 



dip (m 2 a 2 - |) 



[1- 



a cost 



09 



1 



cos 6]' 



■ip+Ptp = Q. (14) 



Eq. (14) is the Schrodinger equation that corresponds to 
He- It is the analog to Eq. (5) wherein the q dependence 
has decoupled from the the angular part of the kinetic en- 
ergy operator and a curvature potential Vc results from 
insisting upon conservation of the norm. 

The Hamiltonian Hq for a particle that lives on the sur- 
face may be obtained by the method employed to derive 
H q by setting q — in Eq. (1) from which 



dr ■ dr — a d9 + (R + a cos6 



(15) 



The resulting expression is simple; Hq is Eq. (11) with Vc 
omitted [2(| . It should be emphasized that for more com- 
plicated surfaces the kinetic energy operator will have 
terms depending on the surface curvature not present 
here because of the azimuthal symmetry of the torus |15| . 
The normalization of an eigenfunction is determined by 



f(lf /°2tt />27r 
Iqi Jo Jo 

with 



tp*(q,9,<P)iP(q,9,(p)M(9,q)d9d(pdq = 1. 

(16) 



M(9,q) =a q F q 



(17) 



when q ^ 0. Wave functions obtained from -f/o.c are 
normalized with M(6, 0) and the q integration omitted. 



ld 2 X , 



2 8q 2 



+ V n (q)Xn = E n Xr, 



with a = a/R,(3 = 2Ea 2 and F = 1 + a cos9. The 
curvature potential Vc is 



V c 



1 1 
ia^F 2 ' 



(12) III. COMPUTATIONAL METHOD 

The goal is to obtain eigenvalues/functions of H q that 
can be compared to those of He and Hq. A procedure for 
determining the low-lyi ng e igenvalues and eigenfunctions 
of Hq has been given in 26] and applied to He in |25[ so 

(13) 

the focus here may be placed on the method employed for 
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solving Eq. (5). In the q — > limit the surface solutions 
are independent of the specific choice of V n (q), but for 
finite q a form or forms for V n (q) must be settled upon. 
Two convenient choices are hard wall confinement with 
the walls at ±L/2 and an oscillator potential V n (q) = 
LU 2 q 2 /2. 

The main complication in solving Eq. (5) ensues from 
the integration measure for the geometry described by 
Eq. (3), which precludes adopting a simple basis set of 
trignometric functions in 9 since they are not in general 
orthogonal over M(0, q). It should also be noted that for 
finite q it is not possible to recover orthogonality by re- 
scaling the basis states because the resulting Hamiltonian 
matrix is not Hermitian. These difficulties can be avoided 
if one performs a two-variable Gram-Schmidt procedure. 

Because the interest here is focused on issues other 
than generating a large number of eigenfunctions, two 
basis functions in q and three in 9 were used. This is a 
reasonable ansatz; the q motion will produce an energy 
spectrum with much larger spacing than the 0, <p motion, 
so two functions in q are sufficient to insure nothin g is 
missed. Additionally, it has been shown previously |26| 
that only a few trigonometric functions in general are 
necessary to accurately describe an eigenfunction on T 2 . 

For hard walls the basis functions with n = (0, 1, 2) are 

$hl = cos(^)cosM) (18) 

*&=sM^)cos(n0) (19) 
and for oscillator confinement 

$°? c = e -^ 2 cos(«0) (20) 

= e-^ 2 1 Hi{yfcq)<x»W). (21) 

For each case the six states computed from the Gram- 
Schmidt procedure are employed to construct the matrix 

= {*T w , osc \H q \K w ,o SC ) (22) 
that yields eigenvalues and wave functions. 



A~ 2 respectively are shown. The dimensionless eigenval- 
ues Pi are found from subtracting the q degree of freedom 
energy (tt 2 /2L 2 or cu/2) from the eigenvalues found from 
the Hamiltonian matrix defined through Eq. (22) and 
multiplying by 2a 2 . The are compared to those found 
in [29] where Vc was included in the T 2 Hamiltonian and 
in |2q where it was not. These results indicate the soft 
constraint quantities are relatively insensitive to differ- 
ing L and w, and are better matched by the spectra of 
|25j . In tables II and III the ground and first excited 
state wave functions for the six cases described above are 
shown. The results illustrate that hard constraint eigen- 
values and eigenfunctions are very good approximations 
to the physically realistic soft constraint values, at least 
for cases where the length scale that determines the sur- 
face energies of the system is near the curvature length 
scale of the device. Here that scale is set by the minor 
radius a; however, in general as the length scale that sets 
local curvature becomes small, Vc increases such that 
matrix elements may become comparable to the 
largest energy in the system. For a disc or strip struc- 
ture the scales can be very different. In the case of a 
disk for example, the energy scale is set by the radius 
of the disk, but a bump or ripples can be placed on the 
disc at much smaller scales [l^, [l3| • Although the results 
here are relatively independent of whether hard wall or 
oscillator confinement was used, it was found that some 
care must be taken with the choice of V n (q). If instead 
of using hard walls at ±L/2, the walls are placed at 
and L, agrement with the hard constraint spectrum is 
lessened, though by only of order ten percent in both the 
eigenvalues and wave function expansion coefficients. A 
possible explanation for this is the sin(n7rg/L) functions 
always vanish on the q — surface so that some terms 
that multiply curvature functions are zero there. 

The basis set expansion employed here comprises two 
functions in the q degree of freedom; for the sake of 
brevity, only angular eigenvalues/eigenfunctions which 
belong to the q ground state wave function have been 
reported. The surface states that correspond to excited 
normal modes lie much higher in energy than the low- 
lying surface excitations dealt with here, but may prove 
important to device modelling as the g-motion becomes 
more diffusive. 



IV. RESULTS 

Toroidal radii R = 500A and a = 250A were chosen on 
the order of structures that have been synthesized 34, 35, 
[2f| , and surface layer widths as set by L, ui within realistic 
values for confinement regions. It should be emphasized 
that the results which follow are very representative; the 
trends exhibited below were found to obtain for larger 
values of R, a, and L as well as for m ^ and negative 
parity states. 

In table I the spectra for hard wall and oscillator con- 
finement potentials with L = 25, 10 A and to = .05, .1 



V. CONCLUSIONS 

The main result of this paper is the good agreement be- 
tween the low-lying spectra and eigenfunctions resulting 
from Hq and those that emerge from H qi and, as im- 
portantly, the relative disagreement that the eigenfunc- 
tions of H q display when contrasted to those of Hq- If 
a two-dimensional approximation is to be adopted for 
curved nanostructures, the results here indicate that for 
cases where the normal excitations are unimportant the 
physics would be better captured with Hq + Vc than with 
H . 
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TABLE I: Ground, first and second excited state eigenvalues 
Pi for the six Hamiltonians relevant to this paper with R = 
500 A and a = 250 A. 





L = 25A 


L = 10A 


uj = .05A" 2 


uj = .lA" 2 


Ref.[25] 


Ref.[26] 


Pa 


-.3405 


-.3406 


-.3489 


-.3488 


-.3511 


.0 


ft 


.6618 


.6610 


.6515 


.6446 


.6386 


1.1223 


P2 


3.7919 


3.7886 


3.7800 


3.7876 


3.6529 


4.0520 



TABLE II: Ground state wave functions; coefficients are nor- 
malized to the constant term in the series to facilitate com- 
parisons. Terms not shown are at least an order of magnitude 
smaller than those given. 







L = 25A 
L = 10A 
uj = .05A~ 2 
a^.lA" 2 
Rcf.[25] 
Ref.[26] 


(1 - .3676 cose + .0693 cos26>)cosff 
(1 - .3675 cosS + .0693 cos26»)cosf§ 
(1 - .3580 cose + .0669 cos26»)e- 025q2 

(1 - .3567 cose + .0654 cos2e)e" 05q2 

1 - .3679 cose + .0784 cos2e 

1 



TABLE III: First excited state wave functions; coefficients 
are normalized to the dominant cose term in the series to 
facilitate comparisons. Terms not shown are at least an order 
of magnitude smaller than those given. 







L = 25A 
L = 10A 
uj = .05A" 2 
a, = .lA- 2 
Rcf.[25] 
Ref.[26] 


(-.0842 + cose - .1369 cos2e)cosff 
(-.0842 + cose - .1370 cos20)cosf3- 
(-.0879 + cose - .1358 cos2e)e- ° 25q2 
(-.0877 + cose - .1362 cos2e)e-' 05q2 
-.0851 + cose - .1540 cos2e 
-.2500 + cose - .0820 cos2e 
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